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ABSTRACT 

Strong encounters between single stars and binaries play a pivotal role in the evolution of star 
clusters. Such encounters can also dramatically modify the orbital parameters of binaries, exchange 
partners in and out of binaries, and are a primary contributor to the rate of physical stellar collisions in 
star clusters. Often, these encounters are studied under the approximation that they happen quickly 
enough and within a small enough volume to be considered isolated from the rest of the cluster. In this 
paper, we study the validity of this assumption through the analysis of a large grid of single - binary 
and binary - binary scattering experiments. For each encounter we evaluate the encounter duration, 
and compare this with the expected time until another single or binary star will join the encounter. 

We find that for lower-mass clusters, similar to typical open clusters in our Galaxy, the percent of 
encounters that will be “interrupted” by an interloping star or binary may be 20-40% (or higher) in the 
core, though for typical globular clusters we expect <1% of encounters to be interrupted. Thus, the 
assumption that strong encounters occur in relative isolation breaks down for certain clusters. Instead, 
many strong encounters develop into more complex “mini-clusters”, which must be accounted for in 
studying, for example, the internal dynamics of star clusters, and the physical stellar collision rate. 

Subject headings: binaries: general — galaxies: star clusters: general — globular clusters: general 
— open clusters and associations: general — stars: kinematics and dynamics — 
methods: numerical 

1. INTRODUCTION 

The evolution of (collisional) star clusters is often con¬ 
ceptualized, at a basic level, as being governed by the 
combination of the long-range cumulative effects of weak 
stellar encounters, known as “two-body relaxation”, and 
the results of short-range strong stellar encounters be¬ 
tween individual stars and binaries. In Monte Carlo 
models for globular cluster (GC) evolution, this assump¬ 
tion is more than a conceptual convenience, and is in¬ 
herent to the functionality of the code (e.g. Spurzem 
& Giersz 1996; Joshi et al. 2000; Vasiliev 2015). Two- 
body relaxation allows stars to gradually exchange en¬ 
ergy, which evolves the cluster towards thermal equilib¬ 
rium, and drives the processes of evaporation, mass seg¬ 
regation and core collapse. It has long been known that 
close encounters with “hard” binaries (i.e., those with 
relatively large binding energy compared to the kinetic 
energies of cluster stars, Heggie 1975) can halt core col¬ 
lapse by donating energy to other stars in the encounter, 
which can be given back to the cluster through two-body 
relaxation processes. This type of strong encounter may 
decrease the semi-major axis of the binary, and indeed, 
strong encounters can modify all orbital parameters of bi¬ 
naries, exchange stars into and out of binaries, and even 
result in physical stellar collisions. Moreover, strong en- 
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counters are a key component to star cluster evolution 
(Hut 1983), and they can alter a binary population from 
its characteristics at birth (e.g. Ivanova et al. 2005; Hur¬ 
ley et al. 2007; Marks et al. 2011; Geller et al. 2013a,b, 
2015; Leigh & Geller 2015). 

As such, the outcomes of single - binary (1+2) and 
binary - binary (2+2) encounters are well studied (e.g. 
Heggie 1975; Hills 1975; Hut & Bahcall 1983; Fregeau 
et al. 2004), and more recently stellar encounters involv¬ 
ing triples are also gaining importance (Leigh & Geller 
2012, 2013). Apart from direct iV-body star cluster sim¬ 
ulations (Aarseth 2003; Wang et al. 2015), it is typical 
to make the simplifying assumption that such encounters 
happen rapidly enough, and within a small enough vol¬ 
ume, that they are effectively isolated from the rest of the 
cluster. With such assumptions, one can run many few- 
body scattering experiments, each involving perhaps 3-6 
stars, to derive statistical cross sections of the outcomes 
of such encounters (e.g. Hut & Bahcall 1983; Fregeau 
et al. 2004), and apply this knowledge to help under¬ 
stand the more complex evolution of a full star cluster, 
which itself may contain many hundreds to millions of 
stars. 

In this paper we investigate the validity of the assump¬ 
tion of treating these strong encounters as isolated. In re¬ 
ality the encounters occur within a star cluster, and most 
often in the dense cluster core, where they may not be 
allowed to progress fully on their own. More specifically, 
we use the numerical scattering code FEWBODY (Fregeau 
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et al. 2004), to perform millions of isolated 1+2 and 2+2 
scattering experiments, described in Section 2. We then 
compare the total encounter duration to the predicted 
time until another single or binary star would join, or 
“interrupt”, this ongoing stellar encounter, as defined in 
Section 3. We find that for certain cluster parameters, 
the fraction of interrupted strong encounters can reach 
over 40% (Section 4). Finally in Sections 5 and 6 we dis¬ 
cuss the implications of these findings and provide our 
conclusions. 

2. Af-BODY SCATTERING SIMULATIONS 

We present results from > 10^ individual 1+2 and 
2+2 numerical scattering experiments performed using 
the FEWBDDY code (Fregeau et al. 2004). We choose pa¬ 
rameters for these scattering experiments that are rele¬ 
vant for the cores of star clusters with total masses (Md) 
and half-mass radii (i?hm) covering the range of observed 
open clusters (OCs) and GCs in our Galaxy. Specifically, 
we sample a grid^ in Md extending from 10^ to 10 ® Mq 
in steps of 0.5 in logiQ^Md [M0] ), and in i?hm from 1 
to 10 pc in steps of 1 pc. We performed two sets of ex¬ 
periments over this grid, one drawing from a mass func¬ 
tion appropriate for an old GC, at an age of 10 Gyr and 
[Fe/H] = —1.5, and the other drawing from a mass func¬ 
tion more appropriate for an OC, at an age of 300 Myr 
and [Fe/H] = 0. 

For a single star in a given scattering experiment, 
we choose a stellar mass from a Kroupa et al. (1993) 
initial mass function (IMF) between 0.1 Mq and the 
turnoff mass at the age of the cluster (^0.95 Mq at 10 
Gyr and [Fe/H] = —1.5, and ^2.96 Mq at 300 Myr 
and [Fe/H] = 0). For a binary, we first choose the 
primary mass (mi) from the same mass function and 
within the same mass limits. We then draw a mass ra¬ 
tio {q = rni/m\) from a uniform distribution to select 
the secondary mass (m2), and enforce the criteria that 
g < 1 and m2 >0.1 M0. We derive a radius for each 
star using this stellar mass and the cluster [Fe/H], follow¬ 
ing the method of Tout et al. (1996), which we provide 
to FEWBODY for identifying physical collisions during the 
encounters (see Fregeau et al. 2004 and Leigh & Geller 
2012 for more details). 

We choose binary orbital elements from the observed 
distributions of binaries with solar-type primary stars in 
the Galactic field from Raghavan et al. (2010), which are 
also consistent with observations of solar-type binaries 
in OCs (e.g., Geller et al. 2010, 2013a; Geller & Mathieu 
2012). Specifically, we draw orbital periods from a log¬ 
normal distribution (with a mean of log(P [days]) = 5.03 
and cr = 2.28), with a short period limit at the Roche 
radius (Eggleton 1983) and a long-period limit at the 
“hard-soft boundary” of the core. Thus we require de¬ 
tached hard binaries initially, and assume that soft bi¬ 
naries are disrupted promptly and are generally not as 
relevant for the dynamical evolution of the cluster. We 
estimate the maximum period for a given hard binary, 
Phs, using the virial theorem, such that, 

ttG /^mim2^^^% , ^-1/2 _3 0^ 

Phs = - mi + m 2 ) cTp , (1) 

V2 \ ma J 

^ As is clear from the top panel of Figure 4, Milky Way star clus¬ 
ters do not occupy this entire parameter grid, though the outlying 
grid points may be useful under other conditions. 



Fig. 1. — Cumulative distributions of the mean sizes (itenc) of 
1-1-2 (solid lines) and 2-1-2 (dashed lines) encounters drawn from 
clusters with ithm = 3pc, over all Md in our grid (as annotated in 
the figure), and using the 10 Gyr mass function. 

where m3 is the mass of the incoming object (either a 
single star or the combined mass of a binary), and cto is 
the three-dimensional velocity dispersion (assumed to be 
(To = V^cto.id). As we discuss in Section 5, we ran addi¬ 
tional experiments drawing binary orbital periods from 
a uniform distribution in log(P), a common theoretical 
assumption. Finally, we draw eccentricities from a uni¬ 
form distribution, and choose all angles of the encounter 
randomly. 

For the majority of our experiments, we choose the im¬ 
pact parameter for a given encounter randomly from a 
uniform distribution between 0 and 1 times the binary 
semi-major axis in a 1+2 encounter, or the sum of the 
two binary semi-major axes for a 2+2 encounter. We dis¬ 
cuss additional experiments with larger impact parame¬ 
ters in Section 5. The velocity at infinity for the incoming 
object is drawn from a lowered Maxwellian distribution 
(typical for star clusters), which is defined by the veloc¬ 
ity dispersion and the escape velocity, both calculated at 
the center of a Plummer (1911) model (with the given 
Mci and i?hm)- 

The above parameters define discrete 1+2 and 2+2 
scattering experiments that are appropriate for the cores 
of typical Milky Way star clusters. For each of our 90 
clusters, we perform 10^ 1+2 and the equivalent number 
of 2+2 unique numerical scattering experiments (at two 
different cluster ages and two different period distribu¬ 
tions, as well as a subset at higher impact parameters). 
We discuss results from these > lO’’ scattering experi¬ 
ments in the following sections. 

3. ENCOUNTER TIMESCALES 

There are two relevant timescales for the purposes of 
this study. The first is the encounter duration, which we 
will call td- The second is the time until another sin¬ 
gle or binary star will join, or “interrupt”, this ongoing 
1+2 or 2+2 encounter, which we will call tg- The first 
timescale, td, is calculated and outputted by FEWBDDY; 
we refer the reader to Fregeau et al. (2004), and specif¬ 
ically Sections 3.3.4 and 3.3.5 and Equations 5 and 6, 
for details. In short, FEWBDDY uses stability assessment 
techniques along with a few simple rules to determine 
when the separate components of the encounter (includ¬ 
ing stable hierarchies) will no longer interact with each 
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other or evolve internally, at which point the encounter is 
deemed complete. These criteria depend on a tolerance 
parameter <5 = Etid/Trei, the ratio of the tidal-to-binding 
(i.e., relative) forces at apocenter, which we set to the 
FEWBDDY default value of 5 = 10“®. (In general, smaller 
values of <5 yield more accurate results.) As we discuss 
in Section 5, our results are only minimally sensitive to 
changes in 6. 

To calculate the second timescale, te, we must first find 
the time for a single (i) or binary ( 2 ) star to interrupt 
an ongoing 1+2 ( 3 ) or 2 - 1-2 ( 4 ) encounter, which we will 
refer to as, respectively, T 3 + 1 , T 3 + 2 , T 4 + 1 , and r 4 + 2 - We 
define these timescales as follows. 

/ lO^ pc-3 \ / \ / 0.5 Mq \ /I AU\ 

V no yV5kms“VV (m) J \ Rene ) 

t-3+1 = 5.4 X 10^° (1 - /b)“^ P yr , 
r 4 +i = 4.2 X 10^° (1 - /b)"^ P yr , 
r3+2 = 4.2 X 10i°/b”'/3 yr , 

T 4+2 = 1.8 X lO^^f^^P yr , 

( 2 ) 

where uq is the central number density, Ui.„is,o is the cen¬ 
tral root-mean-square velocity (and we assume frms.o = 

ao = -v/Sctopd), (w) is the mean mass of a star in the 
cluster, Rene is the mean radius of the encounter, (e.g., 
used to calculate the encounter’s mean geometric cross- 
section, Triigjjj,, see also Figure 1), and /b is the binary 
frequency. We derive these equations following the pro¬ 
cedure in Leonard (1989), with numerical pre-factors de¬ 
rived following Leigh & Sills (2011). Finally, the time 
until a subsequent encounter is given by, 

(bs-i-i + r 3 + 2 ) ^ , for 1 - 1-2 , , 

(r 4 +i -|-r 4 + 2 ) , for 2 - 1 - 2 , 

where T = 1 /r. 

Most of the parameters in Equation 2 come directly 
from the input values discussed above, and the assump¬ 
tion of a Plummer model for the cluster. We describe 
our calculation of the remaining parameters below. 

To calculate no from the Plummer model, we first esti¬ 
mate the total number of stars by randomly drawing suf¬ 
ficient single stars^ from the appropriate mass function 
(discussed above), to reach the cluster mass, Md. In or¬ 
der to reduce stochastic effects, we perform this estimate 
multiple times, until the standard error on the mean 
number of stars derived for the cluster is < O.OlMd. 

We modified FEWBODY to calculate and output the mean 
of the maximum separation between any two stars in the 
encounter at every time step, and divide this value by 
two to get Rene- At a given i?hm, the Rene distribution 
shifts toward lower values with larger Md (see Figure 1). 
This is primarily due to the dependence of the hard-soft 
boundary on Md (through ao), which shifts the distribu¬ 
tions of binary semi-major axes, and also (by definition) 
the impact parameters, towards smaller values at larger 
Md. 

^ This method does not account for binaries, which may not 
have total masses drawn from the IMF, and therefore may slightly 
overestimate the true number density of objects. 


Fig. 2. — Contour plot of the binary frequencies (/p) that result 
from our estimates of the hard-soft boundary, a field binary fre¬ 
quency of 50%, and the 10 Gyr mass functions, for the given i?hm 
and Md values, as discussed in Section 2. (The colors and lines 
show the same values.) 

Finally, to estimate the binary frequency, /b, for a 
cluster with a given Md and i?hm) we first estimate the 
expected average hard-soft boundary in the core of the 
cluster, as the mean of the Phs values (from Equation 1) 
for each encounter. We then calculate the fraction of the 
log-normal input period distribution below this hard-soft 
boundary, fp, and assume that, if the full period distri¬ 
bution could be occupied, this would result in a 50% 
binary frequency (Raghavan et al. 2010). To find the 
(hard) binary frequency for the cluster, we simply multi¬ 
ply /b = 0.5/p. In this way, we assume the total number 
of objects (binaries plus singles) in the cluster is fixed. 
The resulting /b values decrease toward higher Md and 
lower i?hm (Figure 2), as is consistent with observed star 
clusters (e.g., Leigh et al. 2015). 

4. FRACTION OF INTERRUPTED STELLAR 
ENCOUNTERS 

We find from these scattering experiments that the per¬ 
cent of ongoing stellar encounters that are expected to 
be interrupted by an interloping single or binary star in 
a star cluster core ranges from below 1% to over 40%. 
The main results are plotted in Figures 3 and 4. 

For both the 1-1-2 and 2-1-2 experiments, we find that 
te is only minimally sensitive to the cluster mass (Fig¬ 
ure 3). Though the encounter’s geometric cross-section 
(ttR^^P), and binary frequency, both decrease with in¬ 
creasing cluster mass (see Figures 1 and 2), the increase 
in central density roughly cancels out this effect. On the 
other hand, the distribution of to changes dramatically 
with cluster mass, shifting toward shorter encounter du¬ 
rations and developing a pronounced bimodal shape at 
higher Md. 

The shift toward shorter encounter durations at higher 
Md is due to the corresponding decrease in the semi¬ 
major axis at the hard-soft boundary, which results in 
tighter binaries involved in the encounters. The en¬ 
counter duration is then shorter due to the character¬ 
istically lower total angular momentum. 

The emerging bimodal structure in the distribution of 
to at higher Md arises from an increasing frequency of 
physical collisions. Indeed the fraction of encounters that 
result in a physical stellar collision in the cores of our 
most massive and centrally concentrated clusters reaches 
nearly 90% (due to the more compact binaries involved in 
the encounters)^. In the diagonally hatched histograms 

® Note that the collision frequency drops with larger impact pa¬ 
rameters. 
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Fig. 3.— Distributions of the total encounter durations ver¬ 
tically hatched histograms) and times until the next encounter 
(£e; horizontally hatched histograms) for the 1-1-2 (left panel) and 
2-1-2 (right panel) scattering experiments drawn from clusters with 
i?hm = 3pc and the Md values indicated in each panel, for the 10 
Gyr mass function. The diagonally hatched histograms show the 
encounters that resulted in at least 1 (left panel) or 2 (right panel) 
physical stellar collision(s). All histograms are normalized by the 
total number of encounters. 

in Figure 3, we highlight the subset of encounters that 
end with < 2 stars remaining (due to collisions) which 
clearly dominate the peaks at short id- With only two 
(or in some cases one) stars remaining, an encounter is 
considered to be finished, by definition, and therefore 
on average will not last as long as a similar resonant 
encounter that does not result in a physical collision. 

In Figure 4, we compare and tg for each individual 
encounter directly (rather than as an ensemble). In the 
core of a typical GC, at Md ^ 10®'® Mq and i?hm 
3 pc, <1% of encounters are expected to be interrupted. 
However, in the cores of typical OGs, at Md < 10® Mq 
and i?hm < 3 pc, the fraction of ongoing encounters that 
are expected to be interrupted by an incoming single or 
binary star is > 25%, and reaches > 40% for the lowest 
mass and most compact clusters. 

5. DISCUSSION 

We find that a substantial fraction of strong encoun¬ 
ters may be interrupted, while they are ongoing, by an 
interloping single or binary star in the cluster. In this 
case, the common assumption of encounters progressing 
fully in isolation breaks down, and instead the encounter 
develops into a “mini-cluster” with a more complicated 
dynamical evolution. The fraction of interrupted encoun¬ 
ters is largest for OGs (see Figure 4), where, on average, 
the encounter geometric cross-sections are the largest 
(Figure 1) and the encounter durations are the longest 
(Figure 3). For numerical models of star clusters, this 
result is encouraging, as OGs are most often modeled 



2 3 4 5 6 

log,o(M„ [Me]) 


Fig. 4.— Contour plots showing the fraction of interrupted en¬ 
counters (where the encounter duration is greater than the pre¬ 
dicted time until the next encounter: > fe, see Equation 3) 

for clusters with the given total masses, Md, and half-mass radii, 
i?hm- The top panel shows the average fraction from all encoun¬ 
ters, weighted by the 1-1-2 and 2-1-2 encounter timescales from Leigh 
Sc Sills (2011). The bottom and middle panels shows the separate 
results for our 1-1-2 and 2-1-2 scattering experiments, respectively. 
The color shading and solid lines show (the same) results for the 10 
Gyr (and sub-solar metallicity) mass function. In the top panel we 
also show results for the 300 Myr (and solar metallicity) mass func¬ 
tion with the dotted lines. For reference we also plot a sample of ob¬ 
served GCs (Harris 1996) with the gray circles, and OGs (Piskunov 
et al. 2008; Kharchenko et al. 2013) with the black squares. (Note, 
for the OGs, we estimate i?hm from the King core and tidal radii 
from Kharchenko et al. 2013). 


with direct 7V-body methods (e.g., Aarseth 2003), which 
naturally account for stars or binaries interrupting an 
ongoing encounter. GCs are more often modeled with 
Monte Carlo codes (e.g., Chatterjee et al. 2010; Hypki 
& Giersz 2013), where stellar encounters are assumed 
to run to completion in isolation, though direct iV-body 
models are now approaching realistic GCs (Wang et al. 
2015). Fortunately, as Figure 4 indicates, this assump¬ 
tion is much more valid (albeit not perfect) in massive 
GCs. 

In the following, we briefly discuss a few additional pro¬ 
cesses and parameters that may increase or decrease the 
fraction of interrupted encounters, though all of which 
only contribute factors of approximately order unity. 

First, we did not account for any effects of two-body 
relaxation, which leads to mass segregation and the pref¬ 
erential loss of lower-mass single stars from the cluster, 
and thereby tends to increase (m) and /b (e.g., Geller 
et al. 2013a). Increasing (m) decreases tg, and therefore 
tends to increase the fraction of interrupted encounters. 
Increasing /b gives more weight to the 2-1-2 encounters in 
our calculations, and may increase or decrease the overall 
fraction of interrupted encounters accordingly. 

Second, just as stars can be tidally stripped from a 
star cluster by the Galactic potential, stars involved in a 
strong encounter may endure an analogous process ow¬ 
ing to the cluster potential. Following this analogy, we 
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estimate^ an effective tidal radius for each encounter, 

/ M \ 

where Menc is the total stellar mass of the encounter, 
Mcore is the core mass of the star cluster, Ccore is the core 
radius of the cluster, and we assume that the encounter 
occurs at < Iccore from the cluster center. (Here, we es¬ 
timate both Mcore and rcore from the Plummer model.) 
If we assume that any encounter with, i?enc > 
would be tidally disrupted before it is interrupted, this 
decreases the fraction of interrupted encounters by ~15% 
(almost independent of Md and i?hm)- 

We ran additional scattering experiments drawing bi¬ 
nary orbital periods from a uniform distribution in 
log(P), within the same period limits as for the log¬ 
normal period distribution. For the same cluster parame¬ 
ters, this decreases the fraction of interrupted encounters 
by a factor of ~2 over nearly all Md and Rhm, owing to 
the larger fraction of short-period binaries (with smaller 
geometric cross sections) as compared to the more em¬ 
pirically motivated log-normal period distribution. 

Additionally, we ran scattering experiments drawing 
impact parameters from uniform distributions extending 
from zero to 5, 10, 15 and 20 times the size of the binary, 
respectively, and find nearly no change in the fraction of 
interrupted encounters from the results presented above. 

We also ran scattering experiments with 6 values of 
10“^, 10“®, 10“®, 10“'^, and 10“^, respectively. We 
find that the fraction of interrupted encounters decreases 
slightly with smaller 6 (i.e., more accurate outcome clas¬ 
sifications in FEWBDDY), and asymptotes at small 6 toward 
a value of ~ 0.7 times the results at our default 6 = 10“^. 

In closing, one important implication of our results 
may be to increase the expected rate of physical stellar 
collisions between two (or more) stars in star clusters. In 


previous papers, we showed that, for fixed total energy 
and angular momentum, the probability of a physical col¬ 
lision during a stellar encounter increases as the number 
of stars in the encounter increases (Leigh & Geller 2012, 
2015). Thus if additional stars join an ongoing encounter, 
this should increase the probability of a collision. More¬ 
over, collision rates (or cross sections) calculated from 
isolated scattering experiments may only be lower limits 
on the true collision rates in star clusters. 

6. CONCLUSIONS 

Strong stellar encounters in star clusters do not always 
run to completion without another stellar interloper cut¬ 
ting in on the gravitational dance. Indeed, we find that 
~ 20 — 40% of encounters in the cores of OCs may be 
interrupted before completion (see Figure 4). Moreover, 
our results suggest that an assumption of isolated en¬ 
counters is frequently invalid in the OC regime, and in¬ 
stead many encounters may develop into small-iV “mini¬ 
clusters”, though this assumption may still be valid in 
more massive GCs. Finally, the probability of a physi¬ 
cal stellar collision increases with an increasing number 
of stars in an encounter (e.g. Leigh & Geller 2015), and 
therefore these interrupted encounters may enhance the 
rate of stellar collisions in star clusters relative to what 
is assumed from isolated encounters alone. 
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